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Abstract 

How a complex network is connected crucially impacts its dynamics and function 
jl 4|. Percolation, the transition to extensive connectedness upon gradual addition of 
links, was long believed to be continuous Q| but recent numerical evidence on „explo- 



sive percolation" suggests that it might as well be discontinuous if links compete 
for addition. Here we analyze the microscopic mechanisms underlying discontinuous 
percolation processes and reveal a strong impact of single link additions. We show 
that in generic competitive percolation processes, including those displaying explosive 
percolation, single links do not induce a discontinuous gap in the largest cluster size 
in the thermodynamic limit. Nevertheless, our results highlight that for large finite 
systems single links may still induce observable gaps because gap sizes scale weakly 
algebraically with system size. Several essentially macroscopic clusters coexist imme- 
diately before the transition, thus announcing discontinuous percolation. These results 
explain how single links may drastically change macroscopic connectivity in networks 
where links add competitively. 
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Percolation, the transition to large-scale connectedness of networks upon gradual addition 
of links, occurs during growth and evolutionary processes in a large variety of natural, 
technological, and social systems Percolation arises in atomic and molecular solids in 
physics as well as in social, biological and artificial networks 0, Q"[lO|. In the more complex 
of these systems, adding links often is a competitive process. For instance, a human host 
carrying a virus may travel at any given time to one but not to another geographic location 
and therefore infect other people only at one of the places HQ. 

Across all percolating 

systems, once the number of added links exceeds a certain critical value, extensively large 
connected components (clusters) emerge that dominate the system. 

Given the breadth of experimental, numerical, and empirical studies, as well as several 



theoretical results and analytic arguments 



13l-|l6|. percolation was commonly believed to 



exhibit a continuous transition where the relative size of the largest cluster increases con- 
tinuously from zero in the thermodynamic limit once the number of links crosses a certain 
threshold. So recent work by Achlioptas, D'Souza and Spencer came as a surprise because 
it suggested a new class of random percolating systems that exhibit "explosive percolation" 



17l |. Close to some threshold value, the system they considered displays a steep increase 



of the largest cluster size with increasing the number of links; moreover, numerical scaling 
analysis of finite size systems suggests a discontinuous percolation transition. This study 



initiated several follow-up works (e.g. lOl. Il8l-|25|) confirming the original results for a num- 



ber of system modifications. These in particular support that competition in the addition 
of links is crucial; the key mechanisms underlying discontinuous percolation, however, are 
still not well understood and the impact of individual link additions is unknown. 

Gaining one or a few links may have drastic consequences for a network's growth and its 
overall dynamics, depending on whether or not such individual links qualitatively alter the 
global connectivity of a network. For instance, spontaneous activity in developing neural 
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circuits may become persistent after establishing some additional synaptic connections 
2?! . Similarly, during beginning pandemics the specific travel patterns of a single infected 
person may substantially change the number of infecteds on a time scale of months 

Here we identify how microscopic single-link additions impact competitive processes. We 
find that in generic percolation processes, single links do not induce macroscopic gaps in 
the largest cluster size as the system size N ^ oo. Nevertheless, the gap sizes decay weakly 
algebraically as A^"'' with often small /3 such that gaps are essentially macroscopic, i.e. 
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(a) no competition 




(b) edge m=2 




(c) edge m=3 




(d) clique k=3 




Figure 1: Competitive vs. non-competitive percolation processes, (a) Non-competitive Erdos-Renyi 
percolation: new randomly chosen links just add. (b) Edge competition: m = 2 links compete with 
each other and clusters of sizes 4 and 5 win the competition and join to form a new cluster of size 
9 (c) 77T, = 3 links compete with each other. Clusters of sizes 3 and 5 join, (d) Clique competition 
(/c = 3): three links within a clique compete. Clusters of sizes 3 and 5 join. Throughout all panels, 
small disks indicate nodes, solid black lines existing links; large shaded disks indicate clusters 
entering the competition with numbers denoting their sizes; red-dashed lines indicate potentially 
new, competing links; solid red lines indicate actual link added. 



substantially large even for systems of macroscopic size ^ 10^^. Such gaps, induced by 
single links, occur at the point of percolation transitions, are a key signature of discontinuous 
percolation, and are announced by several coexisting, essentially macroscopic clusters. 



The nature of discontinuities in competitive percolation processes 

Consider a family of competitive percolation processes where potentially new links com- 
pete with others for addition (Fig. [1]). Starting with an empty graph of a large number 

of isolated nodes (no links), links sequentially add in competition with others. For edge 
competition, for each single-link addition, m potential links are randomly selected. The link 
for which the sum of the cluster sizes containing their two end-nodes is smallest wins the 
competition and adds. Intra-cluster links are possible; these can only broaden the transition 
compared to disallowing them. For m = 1, this process is non-competitive and identical 
to random Erdos-Renyi percolation jl^, whereas for m = 2 it specializes to the process 
introduced before For all m > 2, this kind of competition promotes that during gradual 
addition of links smaller clusters tend to be connected (to form larger ones) before larger 
clusters grow. With increasing m, the competition becomes more strongly competitive, be- 
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Figure 2: Growth of the largest cluster size Ci as a function of the number of added links L for 
non-competitive [k = 2), competitive, and globally competitive percolation processes for both edge 
(inset) and clique competition (main panel); N = 2^^ nodes, quantities on both axes rescaled by 
system size A^. Single realization is displayed for each percolation process. 

cause more potentially new links actually compete. If m is maximal, all potential links in 
the network compete for addition and we have global competition. 

Taking the sum of cluster sizes in edge competitive processes appears somewhat arbitrary 



because, e.g., taking the product 



28 



or, for that matter, any convex function of the two 



6|. We thus consider also clique competition 



cluster sizes, has similar competitive effects, cf. 
that does not suffer from this ambiguity. For clique competition, randomly draw a fixed 
number k of nodes and connect those two of them contained in the two smallest clusters. 
Here k = 2 describes non-competitive random percolation and for all > 3 competition has 
the same principal effect on changes in cluster sizes as edge competition. We remark that 
for maximal possible k we again have global competition. 

For large finite systems, single realizations of genuinely competitive processes (m > 2 or 
k > 3, cf. Fig. [2]) exhibit macroscopic 0{N) changes in the size Ci of the largest cluster 
Ci. In fact, numerical scaling studies (Supplementary Fig. 1) confirm that the transition 
regime in the plane spanned hj i = L/N and Ci = Ci/N shows an 0{1) change of Ci in a 
region of width Ai that scales as A^"''', 7 > 0, for large (cf. also These results may 
suggest that in the limit of infinite systems there is a discontinuous 0{1) gap in the curve 
characterizing competitive percolation in the i — Ci plane. 
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Further investigating the microscopic dynamics of the transition, however, seeds doubt 
about any such gap. If the largest gap ACmax '■= maXi(Ci(L + 1) — Ci{L)) is macroscopic 
(extensive), 

liinti^>0, (1) 

we call such transitions strongly discontinuous, otherwise weakly discontinuous (see Sup- 
plementary Information for an exact definition). For weakly discontinuous transitions, the 
curve in the i — Ci plane does in fact not exhibit any macroscopic gap in the thermodynamic 
limit. 

Evaluating the largest jump size ACmax from extensive numerical simulations of systems 
up to size 2^^ ^ 6.7 x 10^ already suggests (Fig. [3]) that it scales algebraically as 

=^ ~ (2) 

independent of whether the process is non-competitive, minimally competitive (/c = 3, m = 
2) or exhibits even stronger forms of competition (/c > 4, m > 3). As we find that /3 > 
for all such processes, we have limjv^oo ACmax/A^ = and thus the transitions are all 
weakly discontinuous. The only exception seems to be global competition where we find 
/3 indistinguishable from zero and ACmax/A^ ~ 0.5 > for all (Fig. [1]), suggesting a 
strongly discontinuous transition. The set of all numerical analyses therefore suggests that 
competitive percolation transitions are generically weakly discontinuous, and single links do 
not induce a gap in C\ in the thermodynamic limit — oo. Nevertheless, as the gap sizes 
scale weakly algebraically with system size (|2]) with often small /3 such gaps may still be 
essentially macroscopic, i.e. substantially large even for macroscopic systems of large finite 
size A^. 



The impact of single links 

So how can single links actually impact the dynamics of the transition? For the extreme 
case of global competition, exact analytical arguments reveal the occurrence of macroscopic 
jumps and gives key insights about the nature of transitions in competitive percolation 
processes, that similarly hold for weakly discontinuous transitions (see below): We label all 
existing clusters by Ci and their sizes by Cj = |Ci| where the index i enumerates their size 
rank such that Ci > C2 > . . . > Cy^^^ where z/^ax < A^ denotes the total number of existing 
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Figure 3: Gap sizes ACmax decay algebraically with system size N for weakly discontinuous tran- 
sitions. Green symbols: weakest competition; red symbols: stronger competition, black symbols: 
global competition; blue symbols: no competition {k = 2, m = 1). The symbols indicate average 
values for 50 realizations; error bars indicate 25%-quantiles and reflect system-intrinsic fluctuations. 
Solid grey lines are best least-square fits (slopes, /? = 0.013 {k = 5) /? = 0.021 (ni = 3), f3 = 0.065 
(?Ti = 2) and /3 = 0.095 {k = 3)). Black line shows the analytical curve for global competition, 
where /3 = 0. 

clusters. For global competition each newly added link joins the two smallest clusters in the 
entire system such that C,^^^^ + C;^^^^_i — )■ C. For simplicity of presentation, we choose the 
system size to be a power of 2. This ensures that up to Li = N/2 new links only connect 
1-clusters (isolated nodes) to result in new 2-clusters (two nodes with a single connecting 
link) such that the maximum cluster size stays Ci = 2 for all L < Li. The subsequent N/A 
links each connect 2-clusters to 4-clusters, keeping Ci = 4 until L2 = 3N/4. In general, new 
links added between Ln-i and L„ connect ra-clusters to 2?7,-clusters keeping Ci = 2" where 
Ln = 2" ' ^'^^ ^ — ^'^S2\'^)- the final step, at L = — 1, the remaining two 
y-clusters join and induce the largest gap 



analytically confirming the numerical findings (Figs. [2] and [3]). As a consequence, global 
competition (involving information about the entire system's state for local link addition) 
implies a genuine gap of size 1/2 in the main order parameter Ci. 



AC. 

N 



max 



1 

2' 



(3) 
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For weaker forms of competition, closely related link adding mechanisms control the 
cluster joining dynamics. Inspecting the impact of single link additions onto cluster joining 
dynamics in more detail we identify three distinct mechanisms that may contribute towards 
increasing the size Ci of the current largest cluster in more general competitive processes: 

(i) Largest cluster growth: the largest cluster itself connects to a smaller cluster of size 

Ci < Ci and grows, Ci + Ci — Ci, to stay the largest cluster. 

(ii) overtaking: two smaller clusters of size Ci , Cj < Ci join into one that is larger than 

the current largest cluster, Ci + Cj — ?■ Ci, and the originally largest cluster becomes 
second largest, Ci — t- C2; 

(iii) doubling: if there are several clusters of maximal size Ci = C2 = ■■ = C^ for some 
u > 2, two of these join, Ci + Cj — Ci for some i,j G {1, . . . , u}, creating a new 
largest cluster of twice the size of the original one. 

For each single link addition, we denote the probability for normal cluster growth (i) by 
Pgr. We say that = if the probability of normal cluster growth (i) is zero up to the 
point where only two clusters are left in the system and normal growth is the only remaining 
way the largest cluster could grow at all (see Supplementary Information for a more formal 
definition). 

As we show in the following, an arbitrary percolation process with = necessarily 
exhibits a genuine gap and thus a strongly discontinuous transition, i.e. ACmax/N stays 
positive in the limit of infinitely large system sizes A^. As growth (i) is prohibited, the 
largest cluster size changes either by overtaking (ii) or by doubling (iii). During any such 
percolation process adding a link never more than doubles Ci. As a consequence, there 
is a certain L' such that Ci(L') is larger than A^/3 but not larger than 2N/3. When Ci 
will be overtaken (or doubles) one more time at some L = L' + AL, the cluster previously 
largest becomes the second largest, Ci — C2 (or disappears in case of doubling). Thus it 
is guaranteed that during percolation two clusters of sizes Ci > N/3 and C2 > A^/3 are 
generated which necessarily join at some time L > L'. Therefore, in any such competitive 
process, prohibited growth = implies that the largest gap is macroscopic. 



AC. 

A^ 



max 



1 

> - 



- 3' 



(4) 
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Hence, all competitive percolation processes with = display strongly discontinuous 
transitions with a strong impact of single link additions. As we show in the Supplementary 
Information, such a gap necessarily occurs at or beyond = 1; thus for extremal competition 
with Pgr = the percolation point, where the largest cluster becomes macroscopic, does not 
necessarily coincide with the point where the largest gap occurs. 



Single links induce gaps in large finite systems 

Nevertheless, many weakly discontinuous transitions still exhibit essentially macroscopic 
gaps for large finite systems: We conjecture that competitive percolation processes in nature 
(or engineering or the social world), in particular spatially extended systems with limited 
range interactions shall naturally allow the largest cluster to grow, > (as do all compet- 
itive percolation processes for non-global clique and edge competition) and they generically 



exhibit weakly discontinuous (if not continuous) percolation transitions 29|, |30|. In specific 
limiting models analytic mean field considerations yield 

^ - /3 > 0, (5) 

thus confirming ([2]). For instance, in a model variant where largest cluster joins with the 
smallest available with probability p, and otherwise the two smallest clusters join with 
probability 1 — p we analytically find that (see Supplementary Information for a detailed 
derivation) 

/3 = 1 + (6) 

^ ^log[(l-p)/(2-p)] 21og(2) 

for < p ^ 1 scales roughly linearly with p. 



Notably, if largest cluster growth does not occur, p = 0, we have /3 = and AC^ax/N > 
in the thermodynamic limit, consistent with Eq. ([3]). More importantly, these results show 
that even if the largest cluster may grow the slightest, i.e. for the smallest possible size 
increase with arbitrarily small p > the percolation transition is weakly discontinuous, 
because /3 > as soon as p ^ 0. Direct numerical simulations well agree with our analytical 
prediction ([6]), see Figure HI The finding that /3 > as soon as p > is consistent with 
the above general result that for arbitrarily small probability Pg^ > of cluster growth, the 
percolation transition is already weakly discontinuous, often with small positive exponents f3 
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Figure 4: Weakly discontinuous transition in stochastic mixture of largest cluster growth (with 
probability p) and suppressed growth, (a) Double-logarithmic plot of ACmax/-^ vs. N for different 
p. The slopes of the theoretical mean field prediction (black lines; ranging from /? = 0.036 (p = 0.05) 
to /5 = 0.33 {p = 0.45)) asymptotically well fit the gap sizes obtained by numerical experiments 
(symbols) . (b) Indeed, the theoretically derived exponent /3 ([6|) as a function of p (no fit parameter) 
systematically well predicts those found from fitting the data in (a) (red dots). 



and thus essentially macroscopic gaps in large finite systems (see numerical example below). 
More generally, the results above suggest that any process with non-maximal competition 
{including non-maximal edge competition (m = 2) displaying "explosive percolation" 
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32|) generically displays weakly discontinuous transitions. 



Finite size scaling and coexisting large clusters 

Further extensive numerical scaling analysis reveals that the gaps in the generic com- 
petitive percolation processes we consider indeed occur coincident with the point where the 
largest cluster size is discontinuous (Supplementary Figure 2). Moreover, immediately be- 
fore the transition, not only the largest gap size, but also the second largest cluster, the third 
largest cluster etc. appear essentially macroscopic (Fig. [5]). In particular, the maxmium sec- 
ond largest cluster generically exactly equals the maximum gap size, ACmax = max/, C2{L); 
see Supplemenary Information for a derivation. Thus for small (3 the largest cluster is es- 
sentially non-unique, in contrast to standard continuous percolation transitions. Finally, 



analytical arguments also demonstrate that the percolation strength 
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24| . defined as the 



difference in largest cluster size immediately after and immediately before the gap, exactly 
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Figure 5: Several large clusters coexist at a discontinuous percolation transition, (a) Main panel: 
Simultaneous emergence of the largest cluster of size Ci, the second largest cluster C2, and the 
third largest cluster C3 in a competitive percolation process (clique percolation, k = 4, N = 2^^). 
Inset: blow-up of C2 and C3 in the region around the transition point, (b) The maximal sizes of 
second and third largest clusters as a function of network size indicate that they have the same 
order of magnitude and the same scaling that is moreover identical to that of ACmax- In particular, 
maxL Ci{L) ~ AC^ax ~ N'f^ for i G {2, 3} with /? = 0.095 ± 0.001 for A; = 3 and p = 0.036 ± 0.001 
for k = 5. The maximum gap size ACmax (x) in fact exactly equals the maximum size of the second 
largest cluster. Thus, there is no unique large cluster right at the transition even for very large 
finite systems. 



equals the size of the second largest cluster before the transition, which in turn scales with 
the same exponent /3 as the gap size ([2]). Taken together, single link additions induce sev- 
eral new distinctive features of discontinuous percolation transitions and thus serve as a key 
mechanism controlling competitive percolation processes. 

Interestingly, the so-called /c-cores of the evolving graph, serving as the key example of 



the drastic impact of single links in traditional percolation theory |33l. |34| . exhibit dynamics 
very similar to that for Erdos-Renyi percolation, even for extreme processes with = 0. 
The fc-core of a graph is the largest subgraph with minimum degree at least k. As numerical 
simulations indicate (Supplementary Figure 6), the size of the 2-core increases continuously 
from zero whereas /c-cores for all /c > 3 exhibit a discontinuous jump induced by single 
link additions. These results hold for both Erdos-Renyi as well as competitive percolation 
processes. Even for extreme processes with = the 2-core is still continuous, but with 
the location of the transition moved to larger values compared to the point of percolation. 
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The dynamics of /c-cores is thus very similar for competitive and standard, non-competitive 
percolation processes, in stark contrast to the dynamics of the largest cluster size. This 
is true even though, as shown above, the latter is also strongly influenced by single link 
additions. 

Discussion 

These results explain how microscopic mechanisms of single-link addition control the 
dynamics of the size of the largest cluster and impact the type of transition. In particular, 
the exponent /3 tells in how far single link additions change macroscopic connectivity. For 
generic competitive processes /3 is smaller than for non-competitive ones (see Fig. |3]), but our 
numerical and analytic results indicate that they are still distinct from zero. Only processes 
with global competition or other extreme forms of competition yield /3 = and thus a 
discontinuous gap ACmax induced by single link addition. Others, more generic processes, 
typically exhibit /3 > and thus a weakly discontinuous transition. 

It is important to note that percolation processes with only moderate competition may 
already yield very small positive exponents and thus essentially macroscopic gaps (see Fig.|2]). 
Here we used "essentially macroscopic" to mean that (a) the addition of single links in systems 
of physically large size induces gaps that are of relevant size (substantial fraction of system 
size) and that (b) the gap sizes increase with stronger competition (e.g. increasing k) yielding 
a decreasing exponent P ^ as k ^ N. As a consequence, even processes actually exhibiting 
weakly discontinuous transitions may display large gaps in systems of physically relevant size 
(compare with Fig. [3]). For instance, if /3 = 0.02, a system of macroscopic, but finite size 

= 10^^ exhibits a gap of AC^^JN ~ N'^ ^ 0.35 although formally AC^^JN as 
N oo. For many real processes with already moderate forms of competition, we expect 
exponents /3 close to zero, and thus conjecture that single links may have a strong impact 
onto how such a network becomes connected. 

In summary, our results demonstrate how in competitive percolation, keeping the growth 
rate of the largest cluster small, strengthens the impact of single link additions that merge 
smaller clusters. Growing (i) and overtaking (ii) markedly distinguish the microscopic dy- 
namics in systems exhibiting competitive percolation. The more largest cluster growth is 
suppressed, the more relevant the discontinuous gap becomes in large systems of given finite 
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size. Single link additions may then induce an essentially macroscopic gap even for weakly 
discontinuous transitions if competition is sufficiently strong. 

Interestingly, a protein homology network has recently been identified [l^ that displays 
macroscopic features akin to explosive percolation. Individual links may also induce abrupt 
changes in several other growing networked systems, possibly with severe consequences for 
the systems' dynamics and function (compare to 33l435| ) . For instance, growing one or 
a few additional synaptic connections in a neuronal circuit may strongly alter the global 
connectivity and thus the overall activity of the circuit l^^l; specific infected individuals 
traveling to one but not another location may drastically change the patterns of infectious 
diseases llj; and the macroscopic properties of complex systems exhibiting competitive 



aggregation dynamics of physical or biological units may exhibit abrupt phase transitions 
induced by a small set of specific individual bonds newly established, compare, e.g. 
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371. 



Our study thus does not only provide recipes (by looking for certain competitive cluster 
formation) to identify real systems that could exhibit a (weakly) discontinuous pericolation 
transition, but also shows that and how single link additions in such systems may induce 
discontinuous gaps, and in turn a collective, very abrupt change of structure and dynamics. 

The current study answers how single-link dynamics underlies competitive percolation in 
general, but does not tell how single link additions are actually generated and controlled 
in any given real system. Future work must bridge this gap and establish how competitive 
percolation, and in particular the creation of essentially macroscopic jumps due to single 
link additions, is infiuenced by predefined structure, e.g. for percolation processes on lattices 
and in geometrical or topological confinement occurring in nature 10|, Il8l-|20L |22|. 
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